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Dedicated to Professor Tetsuji Miwa on his sixtieth birthday 

Abstract. We introduce middle convolution for systems of linear differential equa- 
tions with irregular singular points, and we presend a tentative definition of the 
index of rigidity for them. Under some assumption, we show a list of terminal pat- 
terns of irreducible systems of linear differential equations by iterated application 
of middle convolution when the index is positive or zero. 



1. Introduction 

Middle convolution was originally introduced by N. Katz in his book "Rigid local 
systems" |7], and several others studied and reformulated it. Dettweiler and Reiter 
[3l H] defined middle convolution for systems of Fuchsian differential equations written 
as 

Note that Eq. fll.ip has singularities at {ti, . . . ,tr, 00} and they are all regular. On 
the theory of middle convolution, the index of rigidity plays important roles and it is 
preserved by two operation, addition and middle convolution. Addition is related to 
multiplying a function to solutions of differential equations, and middle convolution is 
related to applying Euler's integral transformation to them. Middle convolution may 
change the size of differential equations. It was essentially established by Katz [7J that 
every irreducible system of Fuchsian differential equations whose index of rigidity is 
two is reduced to the system of rank one. The procedure to obtain rank one system 
is called Katz' algorithm. Subsequently it is shown that the system has integral 
representations of solutions which are calculated by following Katz' algorithm. 

In this paper we study middle convolution for systems of linear differential equa- 
tions with irregular singularities, which are written as 

(1.2) 
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If mj = (z 7^ 0) (resp. mo = 0), then the point z = ti (resp. the point z = oo) is 
regular singularity. In particular if rriQ = ■ ■ ■ = niy. = 0, then the system is Fuchsian. 
There are several important linear differential equations which are not Fuchsian, e.g., 
Kummer's confluent hypergeometric equation and Bessel's equation. We want to 
extend the theory of middle convolution to include those equations. 

A pioneering work on middle convolution with irregular singularities was carried 
out by Kawakami [8]. He focused on the differential equations 

(1.3) (^/„-T)^ = Avi>, 

which is called Okubo normal form if T is a diagonal matrix. Eq. fll.3p is called 
generalized Okubo normal form, if T may not be diagonalizable. Kawakami con- 
structed a map from generalized Okubo normal forms to linear differential equations 
with irregular singularities with the condition ttlq = 0, and he showed that the map 
is surjective, i.e. every irreducible equation written as Eq.f ll.2p with the condition 
mo = corresponds to an irreducible equation of generalized Okubo normal form. 
He considered middle convolution by using generalized Okubo normal forms, because 
Euler's integral transformation is on well with (generalized) Okubo normal forms. 
Yamakawa [T3] gave a geometric interpretation to Kawakami's map and investigated 
middle convolution with irregular singularities with the condition mo < 1 by Harnad 
duality. 

In this paper, we construct middle convolution including the case mo 7^ directly, 
i.e. without generalized Okubo normal forms nor Harnad duality. Our construction 
is explicit as we will discuss in sections [2] and |31 

We firstly define convolution matrices, which are compatible with Euler's integral 
transformation (see Theorem 12. ip . The size of convolution matrices is nM, where 

n is the size of given matrices A^j^ in Eq. (11.20 and M = r + Yli=o''^i- '^^^ module 
defined by convolution matrices may not be irreducible. We consider a quotient of 
convolution matrices to obtain irreducible modules in section [3l which leads to the 
definition of middle convolution. We can describe the quotient spaces explicitly, and 
it is an advantage of our construction. We propose a tentative definition of the index 
of rigidity, which is expected to be preserved by middle convolution. 

We study middle convolution further on the case that m^ < 1 and A^*^ is semi-simple 
for all i in section HI In particular, we show that the index of rigidity is preserved 
by middle convolution on this case. By applying middle convolutions and additions 
appropriately, the size of differential equations may be possibly decreased. Under the 
assumption of this section, we show that if the index of rigidity is positive and the 
system of differential equations is irreducible, then the size of differential equations 
can be decreased to one by iterated application of middle convolution and addition. 
Moreover we show a list of terminal patterns obtained by iterated application of 
middle convolution and addition for the case that the index of rigidity is zero. 

We give some comments for future reference in section O 

2. Convolution 

Let A = {Ami, • • • ) ^i^^ • • • , ^0*"^) be a tuple of matrices acting on the finite- 
dimensional vector space V (dimV" = n). The tuple A is attached with the system 
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of differential equations (11. 2p . We denote by (A) tlie algebra generalted by A 
{i = 0,...,r, j = 6ifi, . . .rrii). The (A)-module V (or (A)) is called irreducible if 
there is no proper subspace W{c V) such that A^j^W C W for any 
Set 



(2.1) 



M = r 



nM 



i=0 



i=0 



y{0) ^ Y®mo^ y{i) ^ y(B(mi+l)^ (i > 1) . 



We fix yU G C and define convolution matrices A^j^ {i = 0, . . . ,r, j = Si^, . . .rrii) acting 
on V by 



(2.2) 



U 



u 



(0) 

1 

(1) 



A 



,(0) 



,(1) 



where f m^v ^ G (i' = 0, . . . , r, j' = (5j/_o, • • • ""^j') and u^l ' are given by 



{i') 



(2.3) 



in 



r rn-„ 



i = 0, = j' = j, 



i"=Oi"=5,",o 

i"=Oi"=<5,//,o 

otherwise. 
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Namely 



(2.4) 



A 



(0) 



A. 



(0) 
mo 



V 



(2.5) 



A 
(z > 1) 



A 



(0) 
mo 



>mo-j 



A, 



(0) 



AO) 4(1) 

^1 ^mi 



A, 



(1) 4(2) 



A, 



(r) 



>mo + (mi+l)+---+(mi_i + l) 



'mi- J 

^1 ^0 f''-'" ^mi+l 



\ 



(r) 



where is the unit matrix of size n. We denote the tuple of convolution matrices by 



A = [Amo, ■ • ■ , ^rn'i, . . . ,Aq'). Note that the definition of convolution matrices is 
a straightforward generalization of Dettweiler and Reiter [3] . Then we can show that 
the convolution corresponds to Euler's integral transformation on linear differential 
equations, which is also analogous to Dettweiler and Reiter 



[(0) 7(1) 



Theorem 2.1. Assume that Y 



yn{z) 



is a solution of 
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Lei 7 be a cycle such that -^{r{w)yj{w){z — wY)dw = for any j and any rational 
function r{w). Then the function U defined by 



(2.7) 



(2.8) 



U 



( J^w^ ^yi{w){z — wYdw ^ 
y j^w^''^yn{w){z - wYdw j 



Uf{z) 



/ J^{w-ti) ^ ^yi{w){z - wYdw \ 
y J^{w -tiy^-^yn{w){z - wydw J 



satisfies 



(2.9) 



ITT / mo r rrii a 



i=l i=o 



{z-t,y+^ 



u. 



The proof will be given in our forthcoming paper 

Let 7f is a cycle turning the point w = t anti-clockwise. Then the contour [7^, 7tJ = 
Izltilz^lti" for i = 0,...,r (to = 00) satisfies the condition of Theorem 12. 1[ We 
should also consider cycles which reflect the Stokes phenomena if there exists an 
irregular singularity. 



3. Middle convolution and the index of rigidity 



3.1. Middle convolution. We have defined convolution in section [21 although con- 
volution may not preserve irreducibility. We now consider a quotient of convolution. 
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We define the subspaces /C, and of V — V®^ by 



(3.1) 



(i>l) 



V. 



nii-l 
(i) 



= {0}(C 0°)), ^ = 0^», 



/ ^i^t ^rui-l ■ ■ ■ ^0 \ 

A® ... 4) 

■-. : 

V ••• ^gy 



(i) 







>C'(/^) = <^ 



(0) 
'Jmi 



^^^) = £'(/.), £(0) 
(/i^O) 



i=0 

J^Z ... 

■-. : 

V ••• A 



l(0) 



(0) 





Vo/ 



(0) 



i=0 i=i5i,o 



where = -(A« + . . . + 4"^)). 

Proposition 3.1. H^e /ia?;e if /C C /C, if £(//) C £(//) anc/ if £'(//) C £'(//) /or 
a// i, J. 

Proo/. We show that if £'(//) C £'(/x). Assume that v = {vlnl ■ ■ ■ vfK^n}, . . . v^'^'f e 

£'(/i). Then vf = (i 7^ 0, 0), wj^^ = • • • = = -£ and (w^] . . . vfh)'^ satisfies 
the definition of C'{^). In particular we have 



mo 



(3.2) 



i'=i 



mo 



j'=i 
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Then \ [j ^ 0) is written as 



(3.3) 



if. 



(0) 



(0) 



\ 



mo -J 



(0) 



J 



\ 



J 



and (f, 

J = 1 

,(i) 



(0) 



mo -J 



f^^"* ^ . . . satisfies the definition of C'{fi). Hence Af\ G C'{fi) for 



mo- We also have A^j^v = for i 7^ 0, because Vm^ 



,(^) 



The other cases are shown similarly. 



and 



□ 



We define mCf^iV) to be the (mc^ (A) )-module V®^ /{IC + where mc^(A) is 

the tuple of matrices on V^®^/(/C + £(/i)) whose actions are determined by A, and we 
call it the middle convolution of V with the parameter /x. The following propositions 
are analogues to Dettweiler and Reiter [3], which will be shown in our forthcoming 
paper [12]. 



Proposition 3.2. (i) If fij^O, then /C n C{fi) = {0}. 

If fi = 0, then K, + £'(0) C £(0). 
(Hi) If the (A) -module V is irreducible and fi = 0, then /Cn£'(0) = {0} and dim/C + 
dim£'(0) < n(M- 1). 

Proposition 3.3. Assume that the (A) -module V is irreducible, 
(i) mc^iy) ~ y as {A) -modules. 

(a) The {mCfj,{A)) -module mc^{V) is irreducible and V ~ mc-^{mc^{V)) for any /x. 
3.2. Addition. Let y be a solution of Eq. (ll.2p . Then the function 



r rui 



(3.4) 



/i 



mo (0) 

r' = exp|-v^z^-yy 



1=1 



satisfies the equation 
dY' 

(3.5) 



dz 



i=l j=0 



+ ^JfIn 

{z - ti) 



J+1 



Y'. 



We now define addition for the tuple A = {Amo-, • • • , A^^\ Ami, • • • ; ^0^'*) by 
(3.6) M^(A) = A + /!/„ = {At\ + /x^ J„, . . . , A^K f^^^Q, 

where JI = (/imo, • • • , /if^ /XmL . . . , f^o^) G C*^. 
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3.3. Index of rigidity. Let A = {Ami, • • • v 

acting on V. Set 



A^^^ , Am^ , • • • , Aq' ) be a tuple of matrices 



(3.7) 



and 





... /Iq 


\ 




A« 


... A« 




e End(V^®(™'+^)) 









(z = 0,...,r) 





Ami 


I 







... C-0 


Cm; 


... c« 







^ ■■■ 


Cnii 



We define the index of rigidity by 

r 

(3.9) idx(A) = dim(C«) - (M - 1) (dim(\/))2, 

where M = r + X]i=o"^«- '^'^^ condition = is equivalent to 

k 



0, k = 0, . . . , rrii. 



3=0 



The following proposition is readily obtained by Eq.f l3.10|) : 

Proposition 3.4. The index of rigidity is preserved by addition, i.e. idx(Mjj(A)) = 
idx(A). 

Conjecture 1. If the {A)-module V is irreducible, then the index of rigidity is pre- 
served by middle convolution, i.e. idx(mc^(A)) = idx(A). 

We will prove the conjecture for a special case in section H] (see Proposition 14. 2p . 
We define the local index of rigidity by 

(3.11) idxi(A) = idx,[A« , . . . , 4^] = dim(C«) - {rrn + 1) {dim{V)f. 

Then we have 



(3.12) 



idx(A) = ^idxi(A) + 2(dim(\/))2. 



i=0 



3.4. Example. We consider irreducible systems of differential equations of size two 
written as 



(3.13) 



dY 



(0) , ^0 



(1)" 



Y, Y 



It has an irregular singularity at 2; = oo and a regular singularity at 2; = 0. 

First we consider the case that Af^ is semi-simple. It follows from irreducibility 
that Af'^ is not scalar. By applying addition (i.e. multiplying e'^'^z'^' to the solution 
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Y) and gauge transformation (i.e. multiplying a constant matrix to the solution Y), 
we may assume that A^^ is a diagonal matrix with the eigenvalues and —i^i^ 0) 
and A^f^^ has the eigenvalues and —7. Then we may set 



(3.14) 



(0) ^ f 
-u 



A 



(1) 



—a k 

a(7— o) 



a — 7 



It follows from irreducibility that k ^ and a 7^ 0. By eliminating y2{z) in Eg. (13. 131) . 
we have a second order linear differential equation, 



(3.15) 



avyi 



0. 



If z/ = 1, then Eq. fl3.15p represents the confluent hypergeometric differential equation. 
Eq. fl3.15l) for the case z/ 7^ reduces to the confluent hypergeometric differential 
equation by changing the variable z' = vz. The index of rigidity for Eq. fl3.13l) is two, 
because dim(C^°^) = 4 and dim(C''^^) = 2 which follows from z/ 7^ and A; 7^ 0. 

We investigate middle convolution mc^ for the matrices in Eq. (l3.14p . Convolution 
matrices are given as 



(3.16) 



(0) 



A 



(1) 



/ 



Af^ 





4°^ A'i'+^^h 



-a 



\ 



—V — - a — 7 



\ / 

/ 








k 



\ 



—a + /i 
\ -z/ a -7 + yU / 



The dimension of the space /C(~ /Ci) is one and the space is described as 



(3.17) 



/C 



Ker(A^ 



(l)^ 



C 



/ \ 



k 



The space {fi 7^ 0) is described as 



(3.18) 





10 



KOUICHI TAKEMURA 



Hence the dimension of the space C{a) for the case [i^ a (resp. yU = a) is one (resp. 
two). We concentrate on the case \i = a{j^ 0). A basis of the space C{a) is given by 



(3.19) 



Set 



(3.20) 



Then we have 



(3.21) S-^Af^S 



/ 1 \ 






/ 

a(7 — a) 
ku 

\ 



s 



f 


1 


\ 


1 


a(7 — a) 








ku 


k 


^0 





a J 



( -V 







\ 

-kva 



y 



( « — 7 

1/a 



\ 





y -vjoL a J 



Since the second, the third and the fourth column of the matrix 5* are divisors of the 
quotient space mCaiC^) = (C^)®^/(/C + C{a)), the matrix elements of A^^^ and Af'^ 
appear as (1, l)-elements of Eq.f l3.2ip . Hence Eq.f l3.13p is transformed to 



(3.22) 



dy 

dz 



a 



7 



by the middle convolution mCa- The solutions of Eq. fl3.22p is given hy y = c exp(i^2;)2;"~'^ 
(c: a constant). 

We are going to recover Eq.f l3.13|) from Eq.f l3.22|) and obtain integral representa- 
tions of solutions of Eq. fl3.13p . which arise from the equality mc_oW^Ca = id. We 
apply middle convolution mc-a to Eq. (13.22]) . Then we have 



(3.23) 



dW 

dz 



—V a — 7 








-7 



It follows from Theorem 12.11 that the function 

exp(z/w)w 



(3.24) 



W 



exp(z/w)w"^^^-'^(2; — w)^°'dw 



w 



^dw 



is a solution of Eq.f l3.23p by choosing a cycle C appropriately. For simplicity we 
assume u G M>o- Then we can take cycles C which start from w = —00, move along 
a real axis, turn the point w = z or w = and come back to w = —00. By setting 



(3.25) 



W 



a — 7 

V 



-k 



MIDDLE CONVOLUTION WITH IRREGULAR SINGULARITIES 



11 



we recover the matrices in Eq. f l3.14p such as the function W satisfies Eq. f lS.lSp . Con- 
sequently we obtain integral representations of solutions of Eq.f l3.13|) which are ex- 
pressed as 

(3 26) W= — (^ ^ \( /c^^P('^^)^°~^(^~^)~"^^ ^ 

In particular, the function 

(3.27) y{z)= / exp{uw)w''-^-\z - w^dw 

Jc 

satisfies Eq. fl3.15p . and we obtain integral representations of solutions of the confluent 
hypergeometric differential equation. 

Next we consider the case that Af'^ is nilpotent. Set 



Then it follows from irreducibility that 02,1 7^ 0. The index of rigidity is also two. 
But we cannot reduce to rank one case by applying additions Af'^ Af^ + 0/2, 
A'^^ — 7- Aq^^ + I3I2 and middle convolution mc^, because dim(£'(yu)) < 1 for any 
a,(3,fi, which follows from 02,1 7^ 0. Note that solutions of the differential equations 
determined by Eq. fl3.28p are expressed in terms of Bessel's function by setting z = x^. 



4. The case m,- < 1 for all i 



In this section, we investigate the index of rigidity and middle convolution for the 
case rrii < 1 for all i (see Eq.( 11.2p ). The case = is included to the case = 1 



by setting A^^^ = 0. We assume that A^' is semi-simple for all i and V{= C") is 
irreducible as (A)-module. 

4.1. Index of rigidity. To study the index of rigidity, we investigate Eq. fl3.10p for 
the case rrij = 1. We ignore the superscript Then Eq. fl3.10p is written as 

(4.1) A,C, = C,A^, A^Co - CoA, + AoC^ - C^Ao = 0. 

By the assumption that Ai is semi-simple, we diagonalize Ai as 

\ 



(4.2) 



Write 



p-^AiP 



dll-ni 



























di 7^ dj {i ^ j). 



dklrik / 



(4.3) 


P-'AoP = 


( 


.[1,2] 
4 [2.2] 


■ ^0 

■ ^0 








4[fc,2] 
^0 


.[k,k] 
■ ^0 
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where Aq'"'^ is a x matrix. It follows from AiCi = CiAi that Ci is written as 



(4.4) 



p-^CiP 








... 


\ 







... 














V 




cfi 


/ 



where C| is a ri; x ni matrix. Then 
(4.5) p-\AoC,-C,Ao)P = 



^[2,1] ^[1] _ 



.[l,2]^[2] ^[1].[1,2] 
^[2,2]^[2] 



\ 4[fc.l]^[l] r^W4[fc.l] 4[fc,2]^[2] 
By writing 



^[fc]^[fc,2] 



(4.6) 






r^[1.2] 
Up 

^[2,2] 










--^[fc>2] 





we have 
(4.7) 



P-\A^C,-C^A^)P 


(Ci2 - C?l)Ci'''^ 



/ (rfi - rf2)Co 



[1,2] 



Af^^Cf^ 

^[fc,fc]^[fc] 

\ 



(-r[2]^[2,fc] 



(di - 4)Ci^''^ \ 

"0 



(^2 - 4)4'''^ 



V (4 - t^i)cf (4 - 4)C, 



.[fc,2] 







/ 



CfUi;'^^ and Ci*' 



c/j). Elements of Cq 



It follows from AiC^ - C^Ai + AqCi - CiAq = that A'q'''C[ 

(i^ j) is determined as Cf'^'^ = -(4'''cF^ - CfUf'^Vl^^i - 
are not restricted by relations. Hence the dimension of solutions of Eq. fl4.ip is 

k 

(4.8) Y.^{n,r + dimZ{Al'^^)}, 



1=1 



where Z(A[J''1) = {X e C"'^"' = A^q'^^X}. Let 4,6 be the a x 6 matrix whose 

(i, j)-element is given by Sij, q = (gi, . . . , g^) e (gi H \- Qp = n, qi > ■ ■ ■ > Qp > 

1), A = (Ai, . . . , Ap) G CP. Following Oshima [lO], set 

/ ^ilqi him ■ ■ ■ \ 



[1,1] . 



(4.9) 



^(q; A) 








-^91 ,92 





-^2-^92 


-^92,93 





-^3-^93 
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Every matrix is conjugate to L(q; A) for some q, A. Note that if Aj 7^ Xj {i 7^ j) 
then the matrix L(q; A) is conjugate to the diagonal matrix whose multiphcity of the 



eigenvalue A, is g^. If the matrix Aq''' is conjugate to L{ni^i, . . . ,ni^p^; di^i, 
then the dimension of solutions of Eq. (14. ip is given by 



pi 



(4.10) 



1=1 I i=i 

We denote the type of multiplicities of the matrices {Ai, Aq) which are expressed as 
Eqs.(i2]), <^M and A^'^^ (/ = 1, . . . , fc) is conjugate to L{ni^i, . . .,ni^p^; A/,i, . . . , A/,pJ 

by 

(4.11) (ni, ^2, . . . , Uk) - ((rii,i, . . . , ni,pj, (^2,1, . . . , n2,p2), ■ ■ ■ , (^^fe.i' ■ • ■ > ^Kp^))- 

Note that ni = rii^i + ■ ■ ■ + nl^p^ {I = 1, . . . ,k). Then the local index of rigidity of the 
matrices {Ai,Ao) is calculated as 

k f Pi 

(4.12) 2n^ - E (^')' + 

1=1 I j=i 

If Ai = 0, then k = 1, rii = n and we simplify the notation (rii) — ((ni 1, . . . ,nipj) 
by ( ni^i, . . . ,ni^pj). Note that the notation (^^i,i, . . . ,?^i,pi) was already adapted by 
Kostov [9] and Oshima [lOj for the case of regular singularity. 

By combining Eq. (l4.12p with Eq.( l3.12]) we have the following proposition: 

Proposition 4.1. We assume that itlq = ■ ■ ■ = rrir = 1, A^^^ are semi-simple for 
i = 0, . . . ,r and C" is irreducible as (A) -module. Let 



(4.13) 



(n^*-* > > • • ■ > n^i^ii), n^li > ■ ■ ■ > (j = 1; • • • ) ^'■*'*)) be the type of multiplic- 

3 iP'i 



2 1- ■ ■ 1 ^k(^)) 



ziies 0/ (^1*'', ^o"*)- Then the index of rigidity is equal to 



(4.14) 



idx(A) = EE 

i=0 i=l 



W X2 



— 2rn . 



4.2. Subspace. Next we investigate solutions of the equations 
(4.15) Aivo = 0, Aivi + Aovo = 0, 

for the case that the matrices Ai and Aq are expressed as Eqs. fl4.2p . (14. 3 p and di = 
to understand the subspaces /C^*-* and C{X). Write 



(4.16) 



Vi=P 



-1 



v\ ) 



-1 



\ ^0 / 
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It follows from AiVq = that Vq =0 (i > 2). Hence 

/ 

(4.17) p-\AoVo + A,vi) = 



^0 ^0 \ 



[k] 



,[1] 



G Ker(y4Q '"'^^) and (Z > 2) is determined as —Aq'^^Vq^ /di. The elements of fj^' are 
independent. Hence the dimension of solutions of Eq. fl4.15p is ni + dim(Ker(74Q^'^')). 
If the matrix Aq'^^ is conjugate to L{ni^i, . . . , ni p^; di^i, . . . , and di^i = 0, then 
the dimension of solutions of Eg. (14. 15 p is rii + rii^i. 



,[1] 



4.3. Middle convolution. We now study the matrices for which the middle convo- 
lution is applied. Let i G {1, . . . ,r} and rrii = 1. We investigate the matrices A^^^ 

and for the case that the matrices A^-^^ and A^^ are expressed as Eqs. fl4.2p . (14. 3p . 
di = and C" is irreducible as (A)-module. By changing the order of the direct sum 

V' = (C")®*^, the matrices A^^^ and Aq"* are expressed as 



(4.18) {QP^^'y'A'.^QP' 



M\-i /I (On pffiM 



M\-l A(i)r) TD®M 




p-ly^pffi(M-2) 




_0 

p-l]4pffi(M-2) 





where A = (^Aml ■ ■ ■ A^g, v4m|+i . . . j and Q represents the change of the order of 
the direct sum (C")®*^ Set 



(4.19) 



Af = P 



f 

(d^y'in, 





\ 





P-\ 





\ •■■ (4)-'/n, / 

X = -P-'A* {Ao + filn A) p®(*^-i), 
R = -P-'A*{Ao + fiQP 

Since the j— th row blocks of the matrix P^^{AiAf — In)P zero for j > 2, the 
j-th row blocks of the matrix X' = p-^A^PX + p-^ {Aq + A) P®(*^-i) and 
P~'^AiPR+P~^{Ao+nIn)P are also zero for j > 2. We denote the size of the matrices 



Ai\ Aq^ on the space M 



by n{= nM — dim/C^*^ — dim£(/i)). 



Restrictions of the matrices A^ and A^ to the space IC'^^^ {j ^ i) and £(/i) are 
zero, which follow from the definitions of the spaces. Thus the matrix A^^^ on Ai is 
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diagonalized as 
(4.20) 



/ X 
/ 







M' 



( d2ln2 


V 



I X 
/ 





M' 



dkli 








/ 



where M' 



I X 
/ 



[QP® ) ^A4, / is a unit matrix of the suitable size and the 



dimension of the kernel of ^[^•'lyK is n — n + rii. Since the matrix P ^AiP is diagonal, 
the [/, /]-block of RP^^AiP coincides with that of P~^AiPR, and it is equal to the 
/]-block of -P-\Ao + iiIn)P, if / > 2. Hence 



(4.21) 



I X 
I 



/ X 
/ 



M' 



fx-RP-^A^P nX-XX' 
X' 



p-^AiP 
( Ajf ''1 + 2/x7, 



M' 



"2 



* \ 

* 

* 



where X is expressed as 



+1^1 X 



r[i,i] 



Aq'"^ is the matrix obtained by replacing the domain and the range of ^q'^^ to 
C"^/KerylQ^'^'. It follows from irreducibility that the rank of X is equal to the size of 
Thus, if the matrix is conjugate to 

(4.22) L{nf^ , . . . , nfl ; 0, . . . , 0) © L^n^''^ n<;_^>^ ;-//,..., -/.) 
©L(mi,i,...,mi,p/;Ai,i,...,AiyJ (Aij 0, -/x), 

then we have h — n + ^ > n[ '^^ and the matrix x' is conjugate to 

(4.23) L{n-n + , n^^^ , . . . , n<;_^>^ ; 0, . . . , 0) 

© -L(4°\ . . . , n^pl^;ii, . . . , )U) ® L{mi,i, miy^; Ai,i + /x, . . . , Ai,p; + /x). 



16 
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If the matrix Aq''' (/ > 2) is conjugate to L{mi^i, . . . ,mi^p^; Xi i, . . . , A^^pJ, then the 
matrix Aq''^ + 2^1^ is conjugate to L(m;^i, . . . , mi^p^; Xi^i + 2/x, . . . , A;^p, + 2/i). Hence 



(4.24) idx,(mc^(A)) - idxi(A) = {n - n + m)^ + "^{rij] 



+ (dimZ(Aj,'''^ + + {n-n + + ^(dimZ(4'''] + 2fi)f - 2n^ 



i=2 



i=2 



+ {dunZ(A^,'')y + + ^(dim - 2^^ 

= 2{n — n){2n — rii — rii^i) = 2{n — h){2n — dim/C'-*-'), 

where n\^\ = dim(Kery4g '^'). 

We investigate the matrices Af'^ and A'^^ for the case that tuq = 1, the matrices 
Af^ and Aq"-* are expressed as Eqs. (14.21) . (14. 3 p and di = 0. By changing the order 
of the direct sum V = (C")®^, the matrices Af*^ and A^^'^ + fil are simultaneously 
conjugate to 



(4.25) 



A'; 



(0) 









/i 

4°^ + /xJ ~ ( p-'AiP P-'AoP P-^AP®(^'-^^ 




Aq'' {i 7^ 0) that A^^' and A]^' are simultaneously conjugate to 



where A = ( Am\ ■ ■ ■ A^^^ Ami ■ ■ ■ ) • It follows from similar argument to the case 



i(o) 



|(0) 



(0 



(4.26) 



(0) 



(4.27) 



A 



(0) 







V 



M 



■■■ 0\ 

■•• 

■ ■ • 44fc 

0/ 

* ... * \ 



* A, 



[fc,fc] 



* X J 



where X is expressed as 



X 



^[1,1] 

^0 



III X" 
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Aq^'^' is the matrix obtained by replacing the domain and the range of A^'^^ to 
C"^/Ker(Ao^'^' — fil). If the matrix Aq'^^ is conjugate to 

(4.28) L(nS'^\...,n<f;>^;/i,...,/i)©L(nf\...,4°)^;0,..^ 
©L(mi,i,...,miy^;Ai,i,...,A^p/), (Aij 7^ 0,/i), 

then we have n — n + rif^ > nf^'^ and the matrix X is conjugate to 

(4.29) Lih~n + n[^^ ,nf\..., nf^^ ; -/i, . . . , -/i) 

© L{ni^\. . . , n^'^^^ ; 0, . . . , 0) © L(mi,i, . . . , miy^; Ai,i - /i, . . . , Ai,p/ - /i). 
Hence we also have 

(4.30) idxo(mC/,(A)) - idxo(A) = 2{n - n){2n - dim£(/i)). 

We include the case = to the case mj = 1 by setting A^ ^ = and we have 
n = n(2r + 1) - EI=i dim/C» - dim£(/i). It follows from Eqs.(|lJl]), (ICTD that 

r 

(4.31) idx(mc^(A)) - idx(A) = 2n'^ - 2n^ + ^ {idxi(mc^(A)) - idxi(A)} 

= 2n'^ - 2n'^ + 2{n - n) ^2(r + l)n - ^ dim/C^'^ -dim/:(//)j 

= 2n^ - 2n^ + 2{n ~ n){n + n) = 0. 

Hence the index of rigidity is preserved by application of middle convolution, i.e. 
idx(mc^(A)) = idx(A). Namely we obtain the following proposition. 

Proposition 4.2. If rrii < 1, the matrices A^l^ are semi-simple for all i and (A) is 

irreducible, then the index of rigidity is preserved by application of middle convolution, 
i.e. idx(mc^(A)) = idx(A) for all /x e C. 

4.4. Classification. 

Proposition 4.3. Assume that rrij < 1 (i = 0, . . . , r), A^*^ are semi-simple for all 
i and (A) is irreducible. We identify the case = with the case rrii = 1 and 
= 0. 

(i) //idx(A) = 2, then A is transformed to the rank one matrices by applying addition 
and middle convolution repeatedly. 

(a) // idx(A) = 0, then A is transformed to one of the following cases by applying 
middle convolution and addition repeatedly, where d G Z>i. 

(4.32) Four singularities : {(d, d), {d,d), {d,d), {d,d)}, 
Three singularities : {{d,d,d), {d,d,d), {d,d,d)}, 

{{2d, 2d), {d,d,d,d), {d,d,d,d)}, 
{{3d, 3d), {2d, 2d, 2d), {d,d,d,d,d,d)}, 
{{d,d)~{{d),{d)), {d,d), {d,d)}. 
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Two singularities : {{d, d) — ((rf), (rf)), (d, d) — ((c/), (rf))}, 



{(d, d, d) — ((d), (d), (d)), (d, d, d)}, 

{(d, d, d, d) - ((d), (d), (d), (d)), (2d, 2d)}, 

{(2d, 2d) - ((d, d), (d, d)), (d, d, d, d)}, 

{(3d, 2d) — ((d, d, d), (2d)), (d, d, d, d, d)}, 

{(2d, 2d, 2d) - ((d, d), (d, d), (d, d)), (3d, 3d)}, 

{(3d, 3d, 2d) - ((d, d, d), (d, d, d), (2d)), (4d, 4d)}, 

{(5d, 4d, 3d) - ((d, d, d, d, d), (2d, 2d), (3d)), (6d, 6d)}, 

{(5d, 4d) - ((d, d, d, d, d), (2d, 2d)), (3d, 3d, 3d)}, 

{(3d, 3d) - ((d, d, d), (d, d, d)), (2d, 2d, 2d)}, 

{(5d, 3d) - ((d, d, d, d, d), (3d)), (2d, 2d, 2d, 2d)}, 

{(4d, 3d) - ((2d, 2d), (3d)), (d, d, d, d, d, d, d)}. 



Proof, (i) It is enough to show that the size of matrices can be decreased by appro- 
priate apphcation of addition and middle convolution, because the size of matrices is 
reduced to one by applying many times. 
Let 



(4.33) (<\4^...,nW )- 

{{n[\, . . . \,)), (n2_^i, . . . ,n^\,)), . . . , (?t-[(1) , , • • • > '^[.(l) w ))■, 

(n^*-* > 722'* > ■ ■ ■ > ''^iw) "^j*! ^ ■ ■ ■ ^ ^yP^ multiplicities of {A^l\ Aq-*). 

Note that ri~- \ + ■ ■ ■ + ri~''\i~, = rif' . Then 



(4.34) + J2ir^,r < infy + n« J^ng, = n«(n« + n«). 



There exists a number such that 



(4.35) 




because if not we have contradiction as 



(4.36) 




I ^ j j j' 



By combining with Eg. (14.341) . we have 



(4.37) 



n 



+"!;;,,,) >E('4")^ + EE("S')'- 
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Recall that the index of rigidity is calculated as 
(4.38) idx(A) = f: (Xir^r + E ' 2™'- 

i=0 \ j j j' / 

If ELo ^% + ^ 2rn, then 

(4-39) i2 fE(-?)' + E ^ + ^ 2rn^ 

i=0 \ i i j' / i=0 

which contradicts to idx(A) = 2. Hence 

r 

(4.40) J2{n%+n%^^)>2rn + l. 

We apply addition in order that dim/C*^*-' = n^^^^ + n^^*], for z = 1, . . . , r and the 
dimension of the kernel of A^^ is n[°o). Let fi be the value such that dim£(/i) = n^^o) + 
nj5J) ^. If /i = 0, then it follows from Eq.f OOD that dim/C + dim£'(0) - n{M - 1) = 

X]i=o(^|w + i) ^ ^^'^ ^ ^^'-^ contradicts to Proposition 13.21 (iii). Thus 7^ 0, 
the size of matrices obtained by middle convolution is 

r 

(4.41) (2r + l)n-Eh«+4^i)' 

i=0 

and it is no more than n — 1, which follows from Eq.f l4.40|) . Hence the size can be 
decreased by addition and middle convolution. 

(ii) It is sufficient to consider the case that the size of matrices cannot be decreased 
by addition and middle convolution. Let Z*^*-* be the number which satisfies Eq.f l4.35|) . 

X]i=o(^/W + ^% 1) ^ then the size is decreased by middle convolution (see 
Eq. fl4.4ip ) or the system is reducible (the case /i = 0). Hence Y7i=o(''^% + ^% 1) — 
2rn. If EI=o(^;« < 2™, then we have ELo (S.^f )^ + E,v(ng,)^) < 

2rn^ as Eq.( ]4.39p and it contradicts to idx(A) = 0. Thus Ei=o(^;(0 + 1) ~ ^rn. 
Since idx(A) = 0, all inequalities in Eqs. (14.371) . f l4.39p are equalities. In particular 
we have (nj(]))2 + ^j,{n\ll^.,f = nj(])(n[(]) +n\ll^^), which leads to nj],^^., = nj])^^ for 

1 < j' < pg). Then ((ng)^ + E,(^a,)')n^g, = ^(^g) + = E.ir^f + 

Ej Ej'(^*]')^- Hence if Eq. (14.351) is satisfied, then the inequality in Eq. (l4.35p is re- 
placed by equality. Therefore n((n«)2 + E^.(ng)2) < n«(E,(nf)2 + E, E..,(n«^ 
for all /. It follows from summing up with respect to / that (('T-g)^+Ej(^f])^)^/'^g ~ 
E^(nj*^)^ + ^j'i^^j]')"^ for ^- By setting I = and repeating the discussion 
above, we have 

(4.42) ng=ng = ^, ng(p« + 1) = n« (rf^ + 1), 

Pi 
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for all j G {1, . . . ipf*} and / G {1, . . . , /c*^*^}. Since X]i=o(^i*^ + "^i^i) — 2^^; have 



j=0 



Hence r > 1. 

If n^^*"* = n and p^*'' = 1, then the matrices A^^^ and are scalar, and we may 
omit the singularity corresponding to A^^^ and because they are transformed to 
A^^"^ = Aq^ = by applying addition. If n^^ = n and pf^ > 2, then we have fc^*-* = 1, 
= ■■■ = n^'\,) and 2-(nP+nS^^)/n = 1-l/p?. Thus 1/2 < 2-(nf^+ng)/n < 1 

and the equality holds iff p^^ = 2. Since the matrix A^^ is scalar, this case can be 
regarded as = by applying addition. 

If rii^ 7^ n, then /c*^*^ > 2 and we have n = rif' + ■ ■ ■ + n^,*''^) > k'^^^n^j^^^, n^^l-^ ^ < n^^^^ 
and 

(4.44) 2 ^ ^ = 2 — ^>2 ^>2 p-. 

Hence we have 2— (ra^*^ +n[*^]^)/n > 1 and the equality holds iff /c*^*-* = 2, pj^*'' = P2'* = 1 
and n^*^ = = n/2. 

We consider the case fc'-*-' = 1 for all i. Then ^^^o(-^ ^ ^/Pi^) — 2 and we 
have solution for the cases r = 3 and r = 2. If r = 3, then pf'^ = p^^'^ = 
pT' = p^f^ = 2, i.e. the case {{d,d), {d,d), {d,d), {d,d)} {d = n/2). If r = 2, 
then {Pi'\p^i\p^i^) = (3, 3, 3), (2, 4, 4), (2, 3, 6) or their permutations, i.e. the cases 
{{d, d, d), {d, d, d), {d, d, d)} {d = n/3), {{2d, 2d), {d, d, d, d), {d, d, d, d)} {d = n/4) or 
{{?>d, 3(i), {2d, 2d, 2d), {d, d, d, d, d, d)} {d = n/6). 

We consider the case #{i I A;(*) > 2} > 2. Then it follows from Eq.l^M) that 
#{i I > 2} = 2, fcW = A;(i) = 2, = pf = 1 and = n^^ = n/2 {i = 0,1). 
Hence we obtain the case {{d, d) — {{d), {d)), {d, d) — {{d) , {d))} {d = n/2). 

We consider the case #{i | fcW > 2} = 1 and r > 2. We set A;(°) > 2, A;^^) = 
• ■ ■ = k^"^^ = 1 for simplicity. It follows from Eq.f l4.43i) . 2 — {n^f'' + nf'l)/n > 1 and 
2-{n + n^l)/n > 1/2 {i > 1) that r = 2, nf^ = nf^ = n/2, p^^^ = pf^ = 2 and nf} = 

^1*2 = {i = 1,2). It corresponds to the case {{d,d) — {{d),{d)), {d,d), {d,d)} 
{d'= n/2). 

The remaining case is #{i I A;(^) > 2} = 1 and r = 1. We set > 2, k^^^ = 1 for 
simplicity. By Eq.f l4.43l) . we have 

It follows from pj^^ > 2 and Eq. dCTj) that 2 - (nf^ + riS°J)/n < 3/2 and A;(°) < 4. 

If A;*^*^) = 4, then inequalities in Eq.f l4.44]) must be equalities, pf'^ = p'^^ = pf^ = 
pf^ = 1 and nf'^ = n'^^ = nf^ = nf^ = n/4. We have p^l^ = 2 and n^i^ = n^2^ = n/2. 
Hence we obtain the case {{d,d,d,d) — {{d), {d), {d), {d)), {2d, 2d)} {d = n/4). 
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If = 3, then we have (nf^ + nf})/n < 2/3. Since (nf^ + nf})/n + = 1, 

we have pj^^ = 3 or 2, i.e. (nf^ + n^°l)/n = 2/3 or 1/2. If (nf^ + nf^^)/n = 2/3, 
then pf^ = pf^ = pf^ = 1, nf' = nfl = n/3 {i = 1,2,3), pf^ = 3 and n^^-* = 
n^^ = n^^ = n/3, i.e. the case {{d,d,d) — {{d) , {d) , {d)) , {d,d,d)} {d = n/3). If 
{n^i^ + n^|)/n = 1/2, then we have 

(4.46) 2nS°j(pf^ + 1) = 2ng(pf + 1) = 2nf}{pP + 1) = n 

(0) , (0) , (0) (0) (0) , (0) (0) , (0) (0) 

rSy summing up, we nave ■^[n^ i + ^2 i ^^3 iJ ^1 iPi ' ^2 1P2 ' ^3 1P3 ■ "^mce 

nf''' > n^^ > nf\ we have nfl < n^2i ^ Pi'^ — ^2°^ — ^3°^ pf^ — 2- 

If pf'^ = 2, then p^f^ = p^2^ = 2. If pf'' = 1, then we have 2n^r^l = n^ilip^i '' + 

1) = n^lip'^^ + 1) and 3^3"! = n^ilpf'^ + ^^2*^1^2°^. Hence nfl = nf'l + n^2i 

l/{pf + 1) + l/(pf + 1) = 1/2. Therefore (pf^ + l,pf + l,pf + 1) = (3,3,3), 
(6, 3, 2) or (4, 4, 2). If (pf) + l,pW + i^p(o)^i) ^ 3^ 3)^ ^1^^^ ^(o) ^ ^(o) ^ ^(0) ^ ^^3 

and we obtain the case {{2d, 2d, 2d) — {{d,d), {d,d), {d,d)), {3d, 3d)} {d = n/6). If 
(pS°^ + 1,p5°^ + l,pf^ + 1) = (6,3,2), n!°J = n/12, nf} = n/6 and n^°) = n/4 and 
we obtain the case {{5d, Ad, 3d) — {{d, d, d, d, d), {2d, 2d), {3d)), {Qd, Qd)} {d = n/12). 
If {pf^ + l,pf + 1,pP + 1) = (4,4,2), then nf} = nf} = n/8, nf} = n/4, and we 
obtain the case {{3d, 3d, 2d) — {{d, d, d) , {d, d, d) , {2d)) , {Ad, M) } {d = n/8). 

We investigate the case = 2. It follows from Eq. dCT]) and n = nflpf^ +nflpf^ 
that 

(4.47) nf}{pf + 1) = ng(pf + ^) = - ^) (-S>!°^ + 4>f )• 
By erasing the term p^i^ , we obtain 

(4.48) nS + (p«-l)<i = (rf)-2)4>f. 

Hence p^i^ > 3. It follows from n^f'' > n^2^ that n^i\ < 722°], p^i"* > p^2 \ Pi^^ — 
1 + (p!^^ - l)nf}/nf} = (p!^^ - 2)pf and pf^ < p1^V(p!^^ - 2). We consider the 
case p^i^ = 3. Then pg^'' ^ 3- If ^2°"* = 1? then = and it cannot occur. If 
p'^^ = 2, then n'^l = 2nf'\ and p^f'' = 5. It corresponds to the case {{5d,Ad) — 
{{d,d,d,d,d),{2d,2d)), {3d, 3d, 3d)} {d = n/9). If pf^ = 3, then nf'l = nfl and 
pf^ = 3. It corresponds to the case {{3d, 3d) — {{d, d, d), {d, d,d)), {2d, 2d, 2d)} {d = 
n/6). We consider the case p^i^ > 4. Then pf*^ < p^i'^/{pf^ —2) < 2. If p^^ = 2, 
then p^i '' = 4, nf'l = n^l and pf'^ = 2. It corresponds to the case {{2d, 2d) — 
{{d,d), {d,d)), {d,d,d,d)} {d = n/4). If = 1, then nf} = nf}{pf^ - 3)/{p^^^ - 1) 
and pf'^ = 1 + 4:/{p^^^ — 3). Hence 4 is divisible by p^i '^ — 3 and we have p^^'^ = 
4,5 or 7. If p^l^ = 4, then pf'^ = 5 and 3n^i \ = n'^l. It corresponds to the case 
{{5d, 3d) - {{d, d, d, d, d), {3d)), {2d, 2d, 2d, 2d)} {d = n/8). If p^^^ = 5, then pf^ = 3 
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and 2n[°] = n^2i- corresponds to the case {{3d, 2d) — {{d, d, d), {2d)), {d, d, d, d, d)} 

{d = n/5). If p^^^ = 7, then pf'^ = 2 and 3nf'l = 2n^2 v It corresponds to the case 
{{Ad, 3d) - {{2d, 2d), {3d)), {d, d, d, d, d, d, d)} (d = n/7). 

Thus we have exhausted all the cases. □ 

Remark that some patterns in the list of Proposition 14.31 may not be realized as 
an irreducible system of differential equations. In fact the patterns corresponding to 
Fuchsian systems as 

(4.49) {{d,d), {d,d), {d,d), {d,d)}, 

{{d,d,d), {d,d,d), {d,d,d)}, 
{{2d, 2d), {d,d,d,d), {d,d,d,d)}, 
{{3d, 3d), {2d, 2d, 2d), {d,d,d,d,d,d)}, 

for the case d > 2 (resp. d = 1) cannot be realized (resp. can be realized) as 
irreducible systems, which was established by Kostov [9] and Crawley-Boevey [2]. 



5. Concluding remarks 
We give comments for future reference. 

In this paper we gave a tentative definition of the index of rigidity for systems of 
linear differential equations which may have irregular singularities, and we should 
clarify the correctness (or incorrectness) of our definition. A key point would be 
Conjecture [H which is compatible with Proposition \A.2\ and we should consider the 
case that the coefficient matrices are not semi-simple. 

Crawley-Boevey made a correspondence between systems of Fuchsian differential 
equations and roots of Kac- Moody root systems, which was applied for solving addi- 
tive Deligne-Simpson problem. Boalch pLj gave a generalization of Crawley-Boevey' s 
work to the cases which include an irregular singularity. It would be hopeful to develop 
studies on this direction to understand several properties of differential equations. 

Laplace transformation (or Fourier transformation) has been a powerful tool for 
analysis of differential equations. It is known that Okubo normal form fits well with 
Laplace transformation. In fact, Okubo normal form 

d^ 

(5.1) (x/„-T)— = A^, 



dx 



is transformed to 



(5.2) ^y.fr-^w, 



dz 

which is Birkhoff canonical form of Poincare rank one by Laplace transformation 



(5.3) 



V{z) = / exp{zx)'^{x)dx. 
Jc 
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Assume that the matrices T and A are written as 



(5.4) T 



/ . . . \ / Ali'i] . . . Ali-'^l \ 

^[2,1] ^[2,2] ^[2,fc] 



t2ln2 ••• 

■•. ; 



V • • • / \ Al'^'i] . . . Al'^'^l / 



where ^ tj (i ^ j) and A^*'-'! is a x n^- matrix, and we further assume that A, A'*'*' 
{i = 1, . . . ,k) are semi-simple and Eq. fl5.1l) is irreducible. Then the index of rigidity 
for Okubo normal form (Eq. (l5.ip ) is equal to 

k 

(5.5) (K)' + dim(Z(A[^'^l))) + dim{Z{A)) - n\ 

i=i 

which was described by Haraoka [5] and Yokoyama [H] . On the other hand, the index 
of rigidity of Eq. (15.21) can be calculated as a special case of section H] (see Eq. fl4.14l) ). 
and it is equal to Eq. (15.51) . Hence the index of rigidity in this paper fits well with 
Laplace transformation. 

We now observe an example of Laplace transformation. Set 

(5.6) T=(l A = -f "^'^ + 1 

y y ' y 02,1 «2,2 + 1 

in Eq. (15.21) . which corresponds to Eq. (13.281) . It follows from irreducibility that a2^i 7^ 
0. By (inverse) Laplace transformation, we obtain Eq. (l5.1l) . which is rewritten as 

(r^7\ ^ = /_^/^"2,i ^2,2 + 1 \ _ 1 + 1 ai,2 

^ ' rfx 1 X2 V ) X\ ^2,1 02,2 + 1^/ ' 

and it can be reduced to a scalar differential equation by middle convolution as the 
example in section 13.41 Therefore we should develop a theory of Laplace transforma- 
tion as well as the theory of middle convolution which is based on Euler's integral 
transformation (Theorem 12. ip . 

Several important functions are written as a solution of single differential equations 
of higher order 



(5.8) + ai(^)y("-i) + ■ ■ ■ + an^^{z)y' + a„(z)y = 



where aj(z) (i = 1, . . . ,n) are rational functions which may have poles at prescribed 
points {ti, . . . , tr}- Note that we need to treat delicately on writing Eq. (l5.8p into the 
form of systems of differential equations (II. 2p to reflect the depth of the singulari- 
ties. Oshima [11] formulated middle convolution (Euler's transformation) for single 
differential equations of higher order, and Hiroe [6] studied it for the case that the 
differential equation has an irregular singularity at 2; = 00 and regular singularities. 
Hiroe's result includes a part of the content of section H] in this paper with a differ- 
ent situation. Moreover he clarified a structure of Kac- Moody root system, which is 
based on Boalch's study [T]. Studies on this direction should be developed further. 
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